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ABSTRACT 



Aims. We explore the relation between the charge state of polycyclic aromatic hydrocarbons (PAHs) and the extinction curve mor- 
phology. 

Methods. We fit extinction curves with a dust model including core-mantle spherical particles of mixed chemical composition (sili- 
cate core, sp 2 and sp 3 carbonaceous layers), and an additional molecular component. We use exact methods to calculate the extinction 
due to classical particles and accurate computed absorption spectra of PAHs in different charge states, for the contribution due to the 
molecular component, along five different lines of sight. 

Results. A combination of classical dust particles and mixtures of real PAHs satisfactorily matches the observed interstellar extinction 
curves. Variations of the spectral properties of PAHs in different charge states produce changes consistent with the varying relative 
strengths of the bump and non-linear far-UV rise. 

Key words. (ISM:) dust, extinction — ISM: lines and bands — Ultraviolet: ISM 



1. Introduction 

The nature of interstellar dust has implications for interstellar ex- 
tinction, scattering of starlight, interstellar chemistry, the heating 
of diffuse clouds, the deposition of interstellar ices, and the dy- 
namics of star formation. 

Models of interstellar dust take into account a wide range of 
observational informati on tha t constrains the composition and 
size distribution (Draine 2003). This includes the pattern of ele- 
mental depletions, the variety of observed interstellar extinction 
curves, the distribution of the linear polarisation of starlight, the 
properties of scattered light, absorption and emission features, 
and broadband emission in the visible and infrared regions. All 
models involve a significant fraction of dust in the form of solid 
carbon, some o f which may be in the form of crystalline graphite 
dTielensI 11997). Other forms of carbon have been given vari- 
ous names (amorphous, diamond-like, glassy and hydrogenated 
amorphous carbon, quasi-carbonaceous condensate, soot, yellow 
stuff. . . ) yet these all seem manifestations of carbon in which 
particular valence long-range structures and physical/chemical 
properties (e.g., hydrogen content) are emphasised. 

The strongest interstellar extinction f eature is the br oad ab- 
sorption bump at 217.5 nm discovered bv lStecherl (Il965l) . whose 
carrier is as yet unknown. It has been attributed to a plas- 
mon resonance associated with n — > n * transitions of elec- 
trons in graphite (Stecher & Donn 1965), despite some diffi - 
culties with this interpretation (e. g. Draine & Malh otralll993l) . 
Othe rs have tried to ass ociate the bump with amorphous car- 
bon dMennella et all 1 19981) . although the peak and the width 
of the feature do not seem to be simultaneously well repro- 
duced. It has also been suggested that defect site s at the surface 
of the silicate material could carry the feature dSteel & Dulevl 
1987). The intensity of the bump is known to be unrelated 



to the far-UV rise of the int erstellar extinction curve (ISEC, 
iGreenberg & Chlewickill 19831) . while its FWHM shows an ex- 
tremely lo ose relation to the strength of the nonli near far -UV 
curvature dFitzpatrick & Massal 120071 henceforth |FM07|) . al- 
though with an intrinsic scatter much greater than the observa- 
tional errors. 

Leaving aside specialised hypotheses, almost all the pro- 
posed carriers of the bump seem to require some form of aro- 
matic carbon, either as size-restricted graphite pieces (Draine 
1985) or a s single or stacked Polycyclic Aromatic Hydrocarbon s 
(PAHs) dDulev & Seahrdll998tlLi& Draind200UlDule\l2006l) . 
The possible quanti tative relatio n betw een the bump and PAHs 
was suggested by [Joblin et all d!992l) . who compared labo- 
ratory spectra of mixtures of neutral PAHs and the mean 
galactic ISEC, showing the m to be compatible in spec- 
tral shape. Subsequent papers (Malloci, Mulas & Joblin 2004; 
iMalloci et ail 12005; Ma lloci. Joblin & Mulas 2007al) compared 
theoretical spectra of mixtures of PAHs in different charge states. 

Recent models have a significant fraction of carbon in 
PAH molecule s or clusters, described using approximate aver- 
age properties (Weingartner & Draindl200U iLi & Draiiiefe OOl; 
Zub ko. Dwek. & Arendtl 120041) . In this work we incorporate 
PAHs with their real extinction prope rties in the evolutionary 
model of interstellar dust proposed by I Jones. Dulev & Williami 
( 1990). In the latter model, carbonaceous material is in the form 
of mantles on the surfaces of silicate cores, forming through de- 
position processes conceptually similar to the grain core and icy 
mantle model of dust grains in dark clouds. The core-mantle 
structure and chemical composition of dust particles result from 
the history of the environmental conditions they experienced, 
and responded to, so that the observed extinction along any line 
of sight is an integration of components at different stages of 
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their evolution (ICecchi-Pestellini & Williamsl [19981) . Different 
core-mantle models propose different accretion schemes for the 
carbonace ous mantles, e ither through the cyclic processing of 
ices (as in lLi & Gree nberg 1997) o r through direc t condensation 
(as here, following Jones. Dulev & Williams!! 19901) . Within such 
a framework, PAHs must also respond to the same environmental 
conditions, with chemical changes possibly related to observable 
effects. 

We investigate the relation between the optical prop- 
erties of PAHs in different ionisation states with the 
ISEC. We take advantage of the recent availability of 
state-of-the-art absorption spectra of PAHs in several charge 
states dMalloci. Joblin & Mulasl 2007b), and the extinction 
cross-sections of classical core-mantle dust grains computed 
through an exact approach base d on the multipole expansion 
of the electromagnetic fields (Borgh ese. Denti & Sa iia 2007). 
Synthetic extinction patterns are consistent with the observed 
properties of the carriers of the bump and of the non-linear 
far-UV rise. We propose a quantitative interpretation of the rel- 
ative strengths of these features in terms of varying charge states 
of PAHs. 

2. The model 

While considerable uncertainty remains about the physical na- 
ture of int erstellar dust, in particular its shape and morphol- 
ogy (e. g. IWrighti ri987). spectroscopic evidence and the de- 
pletion data show that silicates and carbons are important 
components of it. Silicates ar e predominantly a morphous but 
with a crystalline component (Li & Draine 2001). Carbons in- 
clude both graphitic, sp 2 , and polymeric, sp 3 , valences (e. g. 
Jones, Duley & Williams] 119901) . They are mostly amorphous, 
with some evidence of the presence of crystalline graphite. 

It is unclear whether the silicates and carbons form distinct 
populations of dust grains or are mix ed within the same particles, 
since both hypotheses fit the ISEC. IChiar et all d2006l) showed 
that the excess polarisation in the 3.4 pm aliphatic C — H stretch 
observed towards the galactic centre is much less than that of 
the 9.7 jt/m silicate feature; they interpreted this to rule out that 
the two features could arise from the same population of dust 
particles. However, this conclusion is not necessarily valid if the 
two features arise from different p arts of the same dust grain 
population, e.g., cores and mantles (Li & Greenberg 2002). 

We use the extinction m odel of Cecchi-Pestellini & Williams 
dl998l) and llati et all d2008l) . i. e. a core + multiple mantle model, 
with dust grains composed of a silicate core and two carbon man- 
tles. The deeper carbonaceous layer represents "older", UV pro- 
cessed sp 2 material, and the surface layer "freshly deposited" 
sp 3 material. For the former we ado pted the BE optical constants 
of lRouleau& Martini dl991l 1 19931) . wh ich are the most repre- 
senta tive for sp 2 amorphous carbon (see lCecchi- Pestellini et alj 
1 19971 for a more detailed discussio n); for sp 3 we used the optical 
consta nts for polymeric carbon of Ashok, Varaprasad, & Birch 
(fl99l . 

Following ICecchi-Pestellini & Williamsl d!998l) and 
llati et al.l d2008l) . the silicate cores are a population of small and 
large spherical grains with a power law size distribution aT q , 
a being the core radius. Silicate grains that are accumulating 
carbon mantle s will all have the same thickness of mantle 
dWhittetl 12002). This component is then completely determined 
by mantle thickness and sp 2 fraction. The contribution of PAHs 
to extinction is represented in the Iati et al. (2008) model by two 
Lorentz profiles, one for the n — > n* transitions, accounting for 
the bump, the second for the cr — > cr* plasmon resonance, whose 
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Fig. 1. Weighted sums of cross-sections of PAHs in different 
charge states, in Mb (10~ 18 cm 2 ) per number of carbon atoms 
Nc- The curves correspond to neutrals (black), anions (dark 
gray), cations (gray), and dications (light gray). The vertical 
dashed lines show the intervals corresponding to the bump and 
the non-linear far-UV rise (see Table[TJ. 



low-energy tail contributes to the far-UV rise. Such a model 
reproduces the observed average ISEC with good accuracy, 
requiring an amount of carbon within the budget av ailable in the 
interstellar medium (for more details see llati et al.ll2008l) . 

For a physical interpretation of the molecular component we 
exploit the database of spectral properties of PAHs computed by 
Malloci, Joblin & Mulas (2007b), including 50 molecules in the 
size range 10-66 C atom s, in the charge state s 0, ±1, and +2. For 
details see lMalloci. Joblin & Mulasl d2007bl) . 

Figure Q] displays the weighted sums of the spectra in our 
sample, weighted with the inverse of the number of C atoms of 
each molecule. The two resonances corresponding to the bump 
and the far-UV rise overlap near 6.5 jum' , and their relative 
intensity varies with charge state, as does the depth of the min- 
imum between them. This is systematic for each molecule and 
depends roughly on charge state and molecular size, yielding 
a 2-dimensional manifold. We show average values for clarity. 
In Table Q] we list, for each sum, the ratio between the inte- 
grated intensities of the bump and of the far-UV rise, and the 
cross-section per C atom at the minimum between them. This 
systematic behaviour hints that PAH charge state is linked to 
morph ology p roperties of the observed ISECs, motivating this 
work: IFM07I showed that the non-linear far-UV rise fluctu- 
ates widely in both strength and slope; our calculations suggest 
that their parametrisation of the ISECs is linked to the relative 
strength and position of the n — » n* and cr — > cr* resonances. 



Table 1. For each sum in Fig. Q] first row: ratio between its 
integrals in the bump (fbump, 2.5 to 6.6 pnT 1 ), and the far-UV 
rise (frfse, 6.6 to 10 jum ) regions; second row: cross-section 
(in Mb per C atom) at the minimum between the n — > n* and 
cr — » cr* resonances. 





Anions 


Neutrals 


Cations 


Dications 


./bump 
/rise 
^~mm 


0.89 
4.18 


0.92 
3.06 


1.01 

2.27 


1.23 
1.55 
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Tabl e 2. In the first part fro m the top, Ay, Rv, C3 (intensity of the bump), and c\ (intensity of the non-linear far-UV rise) from 
Fitzpatrick & Massa| (120071) . Second part: best-fitting parameters of the core + multiple-mantle classical particles model. Their 
meanings are explained in Sect. [2] We also list some physical quantities of astronomical interest derived from them, namely the 
abundances of Si and C locked in this component, relative to total H. Third part: total column density N^ AH of C in PAHs, its 
abundance relative to total H, its fraction in each charge state, and the average charge (in e~ charges) per C atom, as obtained from 
the fit. The reported standard deviations include only statistical errors, not systematic ones d ue to the assumptions in the fit. All 
abundances are computed assuming Nh — 5.9 x 10 21 Ay/Ry, with Nh in cirT 2 and Ay in mag (Whittet 20Q2|). 



Parameter 



HD 

220057 



HD 

197512 



NGC 
457 Pechl3 



HD 

147889 



HD 

93222 



A v 

Rv 

Ca, 



0.56 
2.333 
2.665 
0.374 



0.71 
2.433 
4.136 
0.500 



1.58 
3.106 
2.906 
0.256 



4.44 
4.193 
4.228 
0.525 



N£ AH (10 16 cm- 2 ) 
N£ AH /N H (ppM) 

NPAH^PAH 

NPAH+'^PAH 
NPAH+^PAH 

Charge/N£ AH 



5.92 ± 0.48 

42 
0.62 ± 0.20 
0.34 ±0.18 
0.03 ± 0.06 
0.01 ±0.02 
-0.025 ±0.011 



18.5 ±1.0 

107 
0.76 ± 0.07 
0.15 ±0.05 
0.06 ± 0.03 
0.04 ± 0.03 
-0.043 ± 0.005 



16.4 ± 0.5 

55 
0.75 ± 0.03 
0.25 ± 0.03 
0.00 ±0.01 
0.01 ±0.01 
-0.043 ± 0.002 



108 ±4 
173 
0.75 ±0.10 
0.20 ± 0.09 
0.03 ± 0.02 
0.02 ± 0.03 
-0.041 ± 0.008 



2.17 
5.046 
1.942 
0.250 



w (nm) 


1.25 


1.25 


1.25 


1.25 


1.25 


W sp 2/W 


0.93 


0.94 


0.96 


0.96 


0.98 


(«mto - «max)small ( nm ) 


5-12 


5-12 


5-12 


— 


— 


l^min — ^max /large l^UHJ 


42-250 


42-250 


42-250 


110-250 


110-250 


MRN index 


3.531 


3.531 


3.531 


3.5 


3.5 


N£ US, /N H (ppM)" 


43 


37 


57 


14 


12 


N^ st /N H (ppM) 4 


36 


31 


47 


11 


10 



42.7 ±1.9 

168 
0.65 ± 0.05 
0.16 ±0.05 
0.16 ±0.04 
0.03 ± 0.02 
-0.037 ± 0.004 



N£ Kll /N H (ppM)» 



85 



144 



112 



187 



180 



" Lower limits, since the possible presence of large dust is not included here, see text in Sect. [3] 



3. Results 



We m odelled five extinction curves taken from the sample of 
FM07, with Ry ranging from 2.33 to 5.05, Rv being the ratio of 
total to selective extinction. The fitting procedure provides two 
descriptions for the classical component . For those ISEC s with 
Rv ^ 4, the parameter set exploited by llati et al.l (120081) is ap- 
propriate, while for larger Rv the derived grain distribution is 
significantly different: the small particle component is missing, 
large silicate cores have radii larger than 100 nm, and the index 
q of the power-law size distribution has the standard value of 
3.5. This could indicate depletion of small grains due to coagu- 
lation processes (Stepnik et al. 2003). This would produce a tail 
of very large particles in the size distribution, with a negligible 
effect on the modelled extinction, to which our fit is therefore 
insensitive. The abundances reported in Table [2]for the classical 
dust component do not include the possible contribution of this 
tail of large particles, and thus should be considered with some 
caution as lower limits. For the whole set of ISECs the total man- 
tle thickness is 1 .25 nm, while the sp 2 fraction is almost constant 
between 0.93 and 0.98. This large sp 1 fraction is consistent with 
the short lifetime of the CH bond inferred by laboratory experi- 
ments on interstellar grain analogues (Mennella et al. 2001). The 
presence of a small sp^ fraction resulting from the fit, negligible 
for the extinction in the ISECs considered, is consistent with the 



physical evolutionary model. The full set of model parameters 
obtained from the fit are reported in Table [2] 

To isolate the effect of the molecular component, we sub- 
tract from each modelled ISEC the extinction due to the dust. 
We then fit the residual with a linear combination of computed 
cross-sections of PAHs. We adopt as free parameters the num- 
ber of C atoms in each species, for a total of 200, with the phys- 
ical constraint that they must all be positive. Finally, we sum the 
dust contribution with the best-fitting combination of PAH spec- 
tra to obtain the synthetic extinction curve to be compared with 
observations. As shown in Fig. |2l the fit to the observations is 
quantitatively excellent in all cases, in the whole spectral range: 
the synthetic curves are almost always within the 1 cr region, and 
never further than 2cr from the observations. 

From a geometric point of view, the x 1 hypersurface of this 
fit is a 200-dimensional paraboloid. Its vertex is the formal un- 
constrained solution, its curvature tensor is the inverse of the co- 
variance matrix. The paraboloid is not "bowl shaped", but more 
like a trough, narrow along a small number of directions (< 10) 
and shallow (or even almost flat) in the others. This means that 
some linear combinations of the free parameters are strongly 
constrained by the fit, while others are almost free. 

The most strongly constrained combination of parame- 
ters (i. e. the direction of largest principal curvature of the 
paraboloid) is proportional to the total number of C atoms N^ AH 
in all species, while the combinations that are almost free (i. e. 
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the directions in which the paraboloid is flat) shuffle C atoms 
among species without changing the global charge of the mix- 
ture. This is consistent with the interpretation of the bump and 
far-UV non-linear rise in terms of global properties of the whole 
class of PAHs, rather than a feature associated with some in- 
dividual molecule. To find out to what extent, quantitatively, 
charge, as a global property, is well-determined, we repeated 
the fit 20 times for each ISEC, each time perturbing it by white 
noise, proportionally to the errors reported by FM07. The total 
N^ AH , its fraction in each charge state and average charge were 
computed each time, yielding synthetic statistics stemming from 
the errors in the ISECs. Results are reported in Table [2] of the 
Appendix. The statistical standard deviations obtained show that 
these collective properties of PAHs are indeed well-defined by 
the fit, within the given assumptions. In Fig. [2] we also show 
the extinctions due to core-mantle grains. The curves present 
similar trends: classical dust grains provide an (approximately) 
linear baseline over which the contribution due to the molec- 
ular component lies. This is consistent with previous models 
incorporating PAHs (Li & Draine 200 lj IWeingartner & Draind 
2001: IZubko. Dwek, & A rendt 2004]). Since the contribution of 
the classical component is linear in the spectral region of the 
non-linear far-UV rise, this feature is entirely due to the molec- 
ular component. 

4. Conclusions and future work 

Figure |3]compares the ionisation fractions in Table|2]with the ra- 
tio between the C4 (intensity of the non-linear far-UV rise) and 
C3 (intensity of the bump) parameters of FM07. It is consistent 
with our ansatz that the charge state of PAHs is connected to 
the observed variations in the relative intensity of the two fea- 
tures, even if it is probably not the sole factor determining it. A 
detailed, systematic study will be needed to disentangle the re- 
lations between other specific properties of the PAH population 
and specific parameters of the ISECs. 

This result is not strongly dependent on the model used to 
subtract the dust contribution to extinction: changing the param- 
eters of the classical dust component has effects that are orthog- 
onal to those brought about by the molecular component. For 
example, reasonable changes in the size cut-offs and in the in- 
dex of the power-law size distribution only slightly modify the 
linear baseline on top of which the contribution of the molecular 
component lies (cf. Fig.|2|. 

The validity of the physical interpretation we propose relies 
on the two assumptions that (i) the PAH spectra we used may 
represent the real interstellar molecular component, and (ii) the 
cross-sections do not change in a major way between free-flying 
and clustered PAHs. 

We includ e only fully hydrogenated, unsubstituted PAHs. 
iDulevI J2006I) proposed completely dehydrogenated coronene 
as the carrier of the bump, emphasising that strongly dehydro- 
genated PAHs should display a plasmon resonance near the 
bump position. Hydrogenation is an additional degree of free- 
dom that can respond to environmental conditions and possibly 
play a role, together with charge state, in describing the vari- 
ations observed in the ISEC. Since very little is known about 
the UV spectra of dehydrogenated and substituted PAHs, more 
experimental and theoretical work is needed before this can be 
addressed. 

The PAH spectra we used for the fit are known to slightly 
underestimate the energy of the n — > n* transitions, sim ilarly 
across the whole sample dMalloci. Mulas & Jobl in 2004). This, 
together with the limited number of species used to represent the 
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Fig. 2. Comparison between the synthetic extinction curves, 
with best fitting parameters (continuous line) and the observed 
extinction curves (dotted line FM07), with the 1 cr region in gray 
shade. The dashed and dash-dotted lines are the separate contri- 
butions respectively of the core-mantle grains and of the PAH 
mixtures. 



whole class of interstellar PAHs, produces a systematic error in 
the fit. This is more important for ISECs with a narrow bump, as 
this reduces the number of matching molecules and increases the 
effect of the misplacement of the n — > n* transitions mentioned 
above. Using only half of the molecules at random, the differ- 
ence between the best fit and the ISECs doubles, and absolute 
ionisation fractions change by up to about twice the statistical 
errors displayed in Table [2] These absolute values should there- 
fore be considered with caution, even if some models indeed 
predict a significant population of PAH anions (Tielens 2005) in 
regions including the diffuse galactic ISM. Still, total N^ AH and 
the trends in Fig.[3]are instead almost unaffected by variations in 
the PAH sample used for the fit, since they are due to systematic 
variations of the PAH spectra. Reducing the number of species 
at random worsens the quality of the fits (i.e. the x 1 ) but still 
reproduces the variations in the ISECs considered; on the other 
hand, excluding charge states (e.g. considering only neutrals and 
cations) has a much more dramatic effect, i.e. it becomes impos- 
sible to satisfactorily match the observations. 

As to point (ii), we assumed PAHs to be in free molecu- 
lar fo r m, whereas they could also be clustered dDulev & Seahral 
119981: iRapacioli etafl |2006). Clusters of neutral PAHs are 
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Fig. 3. Comparison of the ratio between the c<\ and ct, parameters 
of IFM07I and the physical parameters derived from our fitting 
procedure for a sample of ISECs. 



weakly bound systems in which the electron densities of indi- 
vidual units (and thus their absorption spectra) are not strongly 
perturbed. The effect of clustering may be more important in 
charged systems, but since this is poorly known it is impossible, 
at present, to consider it in more detail. Our estimation is that 
both clustering and a larger chemical variety would smooth out 
fine structure in the resonances, approaching, in the limit of large 
clusters, bulk optical properties. This would bring the absorption 
spectra of such improved mixtures of PAHs towards better agree- 
ment with observations. This will need to be confirmed by future 
calculations. 

Mixtures of PAHs can accurately account for the UV bump 
and non-linear far-UV rise in extinction, and the PAH charge 
state is linked to the relative intensity of the two features. This 
provides a physical explanation and quantitative relationship of 
the 2-dimensional variations in ISECs described by FMQ7|. In 
this framework, it is therefore not surprising that the bump and 
the non-linear far-UV rise appear to be unrelated from an obser- 
vational point of view: this reflects the variations of the spectral 
properties of PAHs in different charge states. 
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